Exact solutions for vibrational levels of the Morse potential via the 
asymptotic iteration method 
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Abstract 

Exact solutions for vibrational levels of diatomic molecules via the Morse 
potential are obtained by means of the asymptotic iteration method. It is 
shown that, the numerical results for the energy eigenvalues of '^Li2 are all 
in excellent agreement with the ones obtained before. Without any loss of 
generality, other states and molecules could be treated in a similar way. 
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J. Introduction 

Since the appearance of the Schrodinger equation in quantum mechanics, there 
have been continual researches for studying Schrodinger equation with an exactly 
solvable potentials by using differenet methods. The range of potentials for which 
Schrodinger equation can be solved exactly has been extended considerably owing 
to the investigations inspired, for example, by super-symmetric quantum mechanics 
[1], shape invariance [2], the factorization method [3-8], and recently the asymptotic 
iteration method (AIM) [9]. 

In recent years much attention has been focused on AIM. This method reproduces 
exact solutions to many differential equations which are important in applications to 
many problems in physics, such as the equations of Hermite, Laguerre, Legendre, and 
Bessel [9] . The AIM also gives a complete exact solutions of Schrodinger equation for 
Posch- Teller potential, the harmonic oscillator potential, the complex cubic, quartic 
[10], and sextic anharmonic oscillator potentials [11, 12]. Very recently we applied 
the AIM and found the exact eigenvalues for the angular spheroidal wave equation 
[13, 14]. 

Encouraged by its satisfactory performance through comparisons with the other 
methods, we feel tempted to extend AIM to solve exactly the one-dimensional 
Schrodinger equation with the Morse potential [15]. This potential has played an 
important role in many different fields of physics such as molecular physics, solid 
state physic, and chemical physics, etc. This potential has been studied by many 
different approaches such as the standard confluent hypergeometric functions [16], 
the algebraic method [17], the supersymmetric method [18], the coherent states [19], 
the controllability [20], the series solutions with the mass distribution [21], laplace 



transforms [22], etc. 

Fortunately, the one-dimensional Schrodinger equation with the Morse potential 
can be reduced to a second-order homogenous linear differential equation, and we 
therefore directly can make use of AIM to formulate an elegant algebraic approach to 
yield a fairly simple analytic formula which gives rapidly the exact energy eigenval- 
ues with high accuracy. Most importantly, the numerical computation of the Morse 
potential energy eigenvalues using this method is quite simple, fast, and the energy 
eigenvalues were satisfying a simple ordering relation. Therefore, one can unam- 
biguously select the correct starting energy eigenvalue. Moreover, this work shows 
that, the AIM can be a simple alternative approach for computing the vibrational 
levels of the Morse potential. 

In this spirit, this paper is organized as follows. In Sec. 2 the asymptotic 
iteration method for Schrodinger equation with the Morse potential is outlined. 
The analytical expressions for asymptotic iteration method are cast in such a way 
that allows the reader to use them without proceeding into their derivation. In Sec. 
3 we present our numerical results compared with other works, and then we conclude 
and remark therein. 

li Formalism of the asymptotic iteration method for Schrodinger equation 
with the Morse potential 

As an empirical potential, the Morse potential has been one of the most useful and 
convenient model, which gives an excellent qualitative description of the interaction 
between two atoms in a diatomic molecule. As we know the rotational energy of 
a molecule is much smaller than that of its vibrational energy, and therefore, in a 
pure Morse potential model the rotational energy of a molecule has been omitted. 



Hence, the Schrodinger equation for Morse potential with angular momentum i = 
is 



+ V{x) 



^>n{x) = En'^nix), 



2jj,dx'^ 
where En are the energy eigenevalues, and V{x) is the Morse potential function 

V{x) = L)e(e~^^^'^"'^'=^ - 2e"'^(^-^^)). (2) 

Equation (1) is called one-dimensional schrodinger equation, if x is defined on the 
whole line (— oo < a; < +oo), and the eigenf unctions are normalized j_ \'$n{x)\'^dx = 
1. However, for real diatomic molecules x should ranges from to oo. 

De is the dissociation energy, Xe is the equilibrium internuclear distance of a 
diatomic molecules, /i is the reduced mass, and P is an adjustable parameter. Morse 
potential has a minimum value at x = Xe, and it is zero at x = oo. At x = 0, V^(0) 
has a finite value of Dgle^^^" — 26^^") that is positive when /^Xg > ln2. 

Starting with Morse's substitution u = e 2 , we rewrite equation (1) in the 
form 

du^ u du p^h j3^% u^ 

Furthermore, we remove the first derivative by proposing the ansatz 

■^niu) = (^n{u)exp{-p{r)/2y, p'{r) = -. (4) 



Which in turn implies 

dw^ u^ 

where 



1 8;Ui?„ 2 8/iDe 



e„(e„ + l) = ---— ^; 7 = -^^- (6) 
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If we further introduce the frequency 

iuo = PV'^Djfx (7) 

of classical small vibrations about the equilibrium position x = Xe, and express the 
energy parameters in unit huo; that is, 

(e„ + -f 

De = AhuJo; En = Enh^O] ^n = TpJ\2 — ' ^^^ 

it is now more convenient to write the eigenvalue problem in the re-scaled form 

- ^-^^ + ^"(^» + ^) $„(^) + 16AV$„(«) = 32A2<|.„(«). (9) 

In order to guarantee the asymptotic behaviour of this eigenvalue problem when 
u — > oo, and u — > we found that this asymptotic behaviour suggests that 
$„(«) should look like 

$„(«)= «(-+^)e-2^"Vn(«). (10) 

This implies that the function /„(«) will satisfy a second-order homogenous linear 
differential equation of the form 

IjM _ (8A« - ^^^)^M _ (12A + 8A6„ - 32A')Uu) = 0. (11) 

The systematic procedure of the asymptotic iteration method begins now by 
rewriting equation (11) in the following form 

fniu) = Xo{u)f'M + soiu)fn{u), (12) 

where 

Ao(m) = (8 Am- ^^" + ^ )^ and so(u) = (12A + SAe^ - 32A'^). (13) 

u 



The primes of /„(«) in equation (12) denote derivatives with respect to u. 

Now, in order to find a general solution to this equation we rely on the symmetric 
structure of the right hand side of equation (12). Thus, if we differentiate equation 
(12) with respect to u, we obtain 

fniu) = \l{u)f'M + Si{u)fn{u), (14) 

where 

Ai(-u) = Aq(-u) + so(m) + Ao('u), and si{u) = s'q{u) + so{u)\o{u). 
Likewise, the calculations of the second derivative of equation (14) yield 

fn'iu) = Hu)f'M + S2{u)fn{u), (15) 

where 

\2{u) = \[{u) + si{u) + Xo{u)Xi{u), and S2{u) = s[{u) + So{u)Xi{u). 
For {k + I)*'*, and {k + 2)*^ derivatives, k = 1,2, ..., one can obtain 

fi'^'\u) = Xk-i{u)f'M + su-i{u)Uu), (16) 

and 

/(^+2)(«) = Xu{u)f'Su) + su{u)Un), (17) 

respectively, where 

Xk{u) = A^_i(m) + Sk-iiu) + Xo{u)Xk-i{u), and Sk{u) = sl-i(w) + So{u)Xk-iiu).{18) 
The ratio of the {k + 2)*^*, and (A; + 1)*^ derivatives, can be expressed as: 



For sufficiently large k, we can now introduce the "asymptotic" aspect of the method; 
that is, 

Skju) Sk-l{u) _ 
Xk{u) Xk-i{u) 

Thus equation (19) can be reduced to 

^Hfi'-Hu)) = -M^, (21) 

du Xk~i{u) 



which yields 

/f +i)(m) = Cexp (J -^f^du) = C,\u-i{u)exp (J{q{u) + Uu))dv\ , (22) 

where Ci is the integration constant, and the right hand side of equation (22) follows 
from equation (18), and the definition of g. Substituting equation (22) into equation 
(16) we obtain a ffist-order differential equation 



fniu) + Qiu)fniu) = Ciexp ( / {q(u) + \Q(u))du ] , (23) 

which, in turn, yields the general solution to equation (12) 



fniu) = exp - / g{u )du 



ff . ^ , // 



C2 + Ci / exp / {Ao('U ) + 2q{u )}du du 



.(24) 



Here, it should be noted that one can construct the eigenf unctions /„(«) from the 
knowledge of g. 

O Numerical results for the vibrational levels of the Morse potential 

Within the framework of the AIM mentioned in the above section, the energy eigen- 
values of the Morse potential £„ are calculated by means of equation (20). To 
obtain the energy eigenvalues e„, ffist equation (20) is solved for e„ where the 
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iterations should be terminated by imposing a condition Sn{u)= as an approx- 
imation to equation (20). On the other hand, for each iteration, the expression 
^n{u) = Sn{u)Xn-i{u) — s„_i (m) A„(-u) depends on two variables: e„, and u. The 
calculated e„ by means of this condition should, however, be independent of the 
choice of u. Nevertheless, the choice of u is observed to be critical only to speed 
of the convergence to e„, as well as for the stability of the process. In this work 
it is observed that, the best starting value for u is the value at which the effective 
potential of equation (9) takes its minimum value, that is when u = 1. Therefore, 
at the end of the iterations we put u = 1. 

The results of the AIM for e„ with different values of n, yield 



_-3 + 8A _-7 + 8A _-ll + 8A 
^0 — 7: 5 ^1 — , 62 — , [2,0) 



respectively, that means 

-4n - 3 + 8A 



for n = 0,l,2, (26) 



2 

The parameters e„ were calculated by means of 18 iterations only. Therefore, the 
exact energy eigenvalues of the Morse potential En are 

(-2(2n + l) + 8A)2 



64A 



-, for n = 0,l,2, (27) 



For numerical illustration, in table I we calculate the vibrational energies of 
the '^Li2 molecule in the A^T,^ electronic state. The parameters of the respective 
Morse potential are explicitly indicated, with the dissociation energy parameter 
in both units cm~^, and hujo; the later is better suited to compare with the unit 
separation of the corresponding levels of the harmonic oscillator and appreciate 



the decreasing separations in the Morse potential The table includes the energy 
eigenvalues according to Morse's exact solution, E. Ley-Koo et al. [23], and to the 
calculations of this work using AIM. One can also compare the results of this work 
with those of H. Ta§eli [24], and can easily judge the accuracy of the AIM. 

Finally, we would like to emphasize that, within the framework of the AIM, 
we have easily obtained the exact bound state solutions for the one-dimensional 
Schrodinger equation for the Morse potential. 
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Table 1: Energy eigenvalues of Morse potential for '^Li2 in the A^J]^ state with 
De = 8940 cm~^ = 34.997 hcuo, xq = 3.10821, and f3 = 0.616. 



n 



e„ (Morse) [15] 



Sn [23] 



SniAlM) 






-34. 


1 


-33. 


2 


-32. 


3 


-31. 


4 


-30. 


5 


-29. 


6 


-28. 


'~7 

1 


-27. 


8 


-27. 


9 


-26. 


10 


-25. 


11 


-24. 


12 


-23. 


13 


-22. 


14 


-21. 


15 


-21. 


16 


-20. 


17 


-19. 


18 


-18. 


19 


-18. 


20 


-17. 


21 


-16. 


22 


-16. 


23 


-15. 


24 


-14. 



4987869262695313 
5130744534329139 
5416483298483712 
5845091444715536 
6416568973024646 
7130915883411006 
7988132175874689 
8988217850415658 
0131172907033879 
1416997345729385 
2845691166502178 
4417254369352221 
6131686954279552 
7988988921284204 
9989160270366071 
2132201001525260 
4418111114761700 
6846890610075462 
9418539487466475 
2133057746934739 
4990445388480325 
7990702412103161 
1133828817803284 
4419824605580693 
7848689775435389 



-34.4987858673604677 
-33.5130728062414320 
-32.5416466840021528 
-31.5845075006434506 
-30.6416552561655102 
-29.7130899505676886 
-28.7988115838512790 
-27.8988201560141817 
-27.0131156670575088 
-26.1416981169820595 
-25.2845675057870984 
-24.4417238334716096 
-23.6131671000377494 
-22.7988973054838091 
-21.9989144498102718 
-21.2132185330174288 
-20.4418095551052090 
-19.6846875160729802 
-18.9418524159217014 
-18.2133042546505699 
-17.4990430322602215 
-16.7990687487502584 
-16.1133814041206413 
-15.4419809983712994 
-14.7848675315028206 



-34.4987858673600556 
-33.5130728062405367 
-32.5416466840014849 
-31.5845075006429106 
-30.6416552561648174 
-29.7130899505671913 
-28.7988115838500462 
-27.8988201560133717 
-27.0131156670571784 
-26.1416981169814555 
-25.2845675057862103 
-24.4417238334714426 
-23.6131671000371455 
-22.7988973054833259 
-21.9989144498099840 
-21.2132185330171161 
-20.4418095551047223 
-19.6846875160728061 
-18.9418524159213639 
-18.2133042546503994 
-17.4990430322599053 
-16.7990687487498924 
-16.1133814041203536 
-15.4419809983712888 
-14.7848675315027016 
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